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Abstract 

The thermal efficiency of a Sparklet is evaluated by measuring the impulse response of a pendulum subject 
to a single spark discharge. The Sparklet is attached to the end of a pendulum. A laser displacement sensor 
is used to measure the displacement of the pendulum upon discharge. The pendulum motion is a function 
of the fraction of the discharge energy that is channeled into the heating of the gas (i.e., increasing the 
translational-rotational temperature). A theoretical perfect gas model is used to estimate the portion of 
the energy from the heated gas that results in equivalent pendulum displacement as in the experiment. 
The earlier results from multiple runs for different capacitances of C = 3, 5, 10, 20, and 40//F demonstrate 
that the thermal efficiency decreases with higher capacitive discharges. 1 In the current paper, results from 
additional run cases have been included and confirm the previous results. 


1. Introduction 

The Sparklet™, pioneered by The John Hopkins University Applied 
Physics laboratory, is an electromagnetic local flow/flight control 
(ELFC) device. 2-8 The basic concept is illustrated in Fig. 1 . A spark 
is discharged within a typical volume of several cm 3 (Stage 1). The 
high pressure gas discharges through a converging nozzle thereby 
generating an impulse (Stage 2). Provided there is a mechanism for 
recharging the gas in the cavity (Stage 3), the sequence can be re- 
peated. Research on Sparklet has also been conducted at the Univer- 
sity of Texas at Austin, 9 ONERA 10 and Rutgers University. 11 

The objective of the paper is to determine the thermal energy 
efficiency of a particular design of a Sparklet. The thermal energy 
efficiency is the fraction of the electrical energy which results in ef- 
fective heating of the gas in the cavity. 

2. Experiment 

A brief description of the experimental apparatus is provided in the following sections. All the details are presented in 
the earlier paper by Golbabaei-Asl et al. 1 

2.1 Pendulum 

To measure the impulse response of a single spark discharge, a prototype pendulum apparatus and Sparklet actuator 
was constructed. 1 The schematics of the pendulum apparatus is shown in Fig. 2 and the components include A. screw 
terminal block, B. connecting chains and wires, C. sharp carbide pins, D. stationary glass bearing surface, E. wood 
dowel, F. Sparklet actuator. The pendulum was a 12.7 mm diameter wooden (Maple) dowel mounted to a pivot block 
assembly on one end and to the Sparklet actuator on the opposite end. The pivot assembly consisted of a rectangular 
block of UHMW polyethylene (19 x 28 x 63 mm) with two sharp carbide scriber pins (from Moody Tools Model 
51-1520) resting on hard glass surfaces. 



Stage 1: Stage 2: Stage 3: 

Energy Discharge Recovery 

Deposition 


Figure 1: Sparklet 4 
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The mass moment of inertia for a compound pendulum can be 
computed from the total rotating mass m (190.43 gram), the distance 
from the pivot axis to the center of the mass L cg (291.4 mm), and 
the pendulum period T (1.5 s), i.e., I = mgL cg ( ^) 2 kg-m 2 . A low 
torque electrical coupling was fabricated by suspending low stiffness 
conductors between two screw- terminal blocks, one attached to the 
stationary wooden support stand and the other to the rotating pivot 
block (See Fig. 2). 



2.2 Sparkjet Actuator 


Figure 2: Pendulum Apparatus 


A capacitive-discharge, ionization-triggered, tri-electrode 
Sparklet actuator was constructed . 1 Fig. 3 presents the 
schematics of the electrical circuit connected to the Sparklet. 
An ionization trigger was used to discharge the capacitor by in- 
troducing an approximately 1 microsecond pulse, V p , into the 
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step up transformer T to create a 25 kV pulse between nodes 5 
and 3. The pulse was provided by a low cost, commercial trig- 
ger source (Information Unlimited, Model Ignitor 10). Resistor 
R 2 helps to isolate the transformer from the spark gap while 
still allowing the ionization pulse to function. The circuit will 
function without R 2 but the influence of the transformer will 
be evident in the loop voltage measurements, particularly V \2 
and V 34 complicating interpretation of the data. Diode D pre- 
vents reverse current and a second, echo discharge due to the 
collapsing magnetic field. 

The high voltage trigger pulse creates an ionization 

channel between nodes 5 and 3 that subsequently allows the capacitor to discharge along path 1 - 2 - 5 - 3 (Fig. 4(b)) 
Fig. 4(c) shows the construction of the actuator . 1 


Figure 3: Sparklet electrical 

Vdc - 550 volt unregulated, V p 


circuit schematic, 
trigger pulse from 
pulse generator, T pulse step up transformer, R\ = 
5 kOhm, R 2 = 1 kOhm, D diode (IN4007-30), C ca- 
pacitor (3, 5, 10, 20, 40yuF) 
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(a) Ionization trigger 



(b) Spark discharge 
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(c) Dimensioned cross-sectional diagram of Sparklet. All 
dimensions in millimeters. 


Figure 4: Sparklet: a,b) photograph and c) diagram 


2.3 Diagnostics 

Since the maximum expected displacement of the pendulum was only a few tens of micrometers, it is important to 
accurately measure the displacement time history. A laser-based optical displacement sensor, Microepsilon optoNCDT 
model ILD-2200, was selected based on its range (2 mm), resolution (0.03 nm), stand-off distance (24 mm) and 
frequency response (20 kHz at -3 dB). The sensor has both analog and digital signal outputs. The analog output was 
used in this study . 1 

The instantaneous voltage across the gap and instantaneous loop current are measured in order to determine the 
input energy required to the main discharge spark. A large voltage drop is expected across both the supply and return 
legs of the wiring since the instantaneous current can be very high (on the order of 1000 amp). It was observed for a 
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particular implementation that roughly two-thirds of the capacitor voltage is dropped across the supply wiring. For this 
reason, a separate pair of voltage monitoring leads was attached directly to the actuator electrodes close to the spark 

gap- 1 


3. Theoretical Model 


3.1 Model Description 


The geometry of the pendulum is shown in Fig. 5. Prior to discharge and as a result of ambient noise the pendulum 
oscillates with small amplitudes. The discharge creates an air jet which generates a force and moves the pendulum. The 
geometry of the cavity is shown in Figure 6. The cavity is initially at ambient pressure poo and ambient temperature 
T^. The energy deposited to the cavity, Q, is assumed to increase the temperature to T to and pressure to p to at a 
constant density (p to = poo) instantaneously at time to. An inviscid perfect gas model is assumed for the air flow to 
avoid elaborate computations required to account for thermochemical reactions. 
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Figure 5: Pendulum configuration 





3.2 Governing Equations 


3.2.1 Pendulum Motion 

Assuming that 1) the pendulum motion is planar, 2) the dimension of the Sparklet is small compared to the pendulum 
arm R, 3) the mass of the gas inside the Sparklet is negligible compared to the mass of the Sparklet (i.e., m s (t) is 
constant), and 4) 6 <$: 1 radian; the pendulum motion can be simplified from Newton’s second law 


where 


d 2 0 
dt 2 


+ oj z 6 = 
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( 2 ) 


where m p is the mass of the pendulum (assuming a uniform density and constant cross section), m s (t ) is the mass of the 
Sparklet, R is the distance from the pendulum pivot point to the cavity center, F(t) is the force from the cavity exhaust, 
and I is the moment of inertia of the pendulum plus cavity system, i.e.. 


r R 

I - I r 2 dm p + m s (t)R 2 (3) 
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3.2.2 Cavity Jet Flow 

The initial total temperature of the gas in the cavity is 


T to — + 


dQ 

Poo F c v 


(4) 
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where 77 g is the effective heat added to the gas resulting in the increase in the translational-rotational temperature 
of the gas), is the initial density of the gas in the cavity, V is the cavity volume and c v is the specific heat at constant 
volume (see Fig. 6 ). The the subscript t Q indicates the total property of the gas at the time to. The quantity 77 is denoted 
the thermal efficiency of the Sparklet. The thermal efficiency in percent is 100 X 77 . The above expression may be 
rewritten as 

T to = Toofl + (y - 1)6] (5) 


where 

e = _nQ_ 

PooV 


(6) 


For a physical interpretation of 6, see Anderson and Knight. 12 

Anderson and Knight have developed an analytical model to predict F{t) in Eq. (1) for inviscid and perfect gas. 12 
The force from the cavity is given by Eq. (7) assuming that the flow is initially sonic at the throat (see Anderson and 
Knight 12 ). 



0 

for 

t < to 


Fit) 

c(r)[i + (y- i)^][i +Pit-t 0 )]~ : ^ 1 - 1 

for 

to < t < to + A h 

(7) 

PooA 


for 

t > to + A t\ 


where Q is the energy deposited into the cavity, V is the cavity volume, is the ambient pressure, A is the cavity 
exit area, and 77 is the thermal efficiency. The energy deposition is assumed to occur instantaneously. The constants are 
given by the following relations 
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and the initial temperature in the cavity can be found from Eq. (5) 

The condition for initially sonic flow at the cavity exit (see Anderson and Knight 12 ) is 
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* 2.23 for y = 1.4 


The time interval A t\ of sonic flow at the cavity exit is 
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The total pressure ratio to ambient pressure during subsonic flow at the cavity exit is determined by the following 
implicit equation 12 


„ (2/ + 5) + 3 log („ + - 3 946 + (4) (+) +++ 


= 0 for y = 1.4 (14) 


where 
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For initially subsonic throat, e < 2.23, (see Anderson and Knight 12 ) 
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The total pressure ratio to ambient pressure for subsonic throat is found by the following implicit relation 
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where p is found by Eq. (15) and 
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The total discharge time for initially sonic throat, e > 2.23, is 


t f =Aj [8.99(1 + 0.4e) 1/7 - 8.64] 


( 20 ) 


and the total discharge time for initially subsonic throat, e < 2.23, is 
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where e and y- are given by Eq. (5) and (19), respectively. 


3.3 Numerical Method 

The energy efficiency r\ is determined by choosing the value that creates a match between theoretical displacement and 
experimental displacement. The governing equation (Eq. (1)) is second order in time and is solved by a fourth order 
Runge-Kutta scheme. 13 The scheme incorporates variable time steps for different phases of pre-discharge, discharge 
(sonic and subsonic throat intervals), and post-discharge. Newton’s method is used to find the numerical roots of the 
implicit relations (Eq. (14) and (17)) to determine F(t) . The two initial conditions, i.e., the angle 6 and the angular 
velocity dO/dt at the instant of spark discharge are obtained from the experiment as well as the natural frequency oj. 


4 . Results 

Experimental and theoretical results for one run of each capacitive discharges, i.e., C = 3,5, 10,20, and 40 pF are 
shown in Fig. 7. The graphs present the pendulum displacement upon a single Sparklet discharge. The values of 
computed efficiency is given for all the cases (Fig. 7). It is worth to note that the experimental and computed results 
are particular of this specific Sparklet design. 

The average efficiency and standard deviation values of multiple runs (5 to 10) for different capacitances are 
summarized in Fig. 8(a). The average efficiency and standard deviation values for different values of dimensionless 
deposited energy are plotted in Fig. 8(b). The results indicate lower efficiencies with higher capacitances, i.e., capacitive 
discharges which lead to initially sonic throat conditions in the cavity. 
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(a) 7j = 8.91% and Q = 181 mJ 


(b) rj = 7.82% and Q = 285 mJ 




(c) 77 = 5.64% and Q = 572 mJ 


(d) 77 = 4.72% and Q = 1130 mJ 



(e) 77 = 3.97% and Q = 1901 mJ 
Figure 7 : Experiment vs theory 
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(a) 77 and e vs. C 



(b) 77 vs. e 


Figure 8: Thermal efficiency and dimensionless energy deposition parameter 

5. Conclusion 

A novel method for determining the thermal efficiency of the Sparklet is proposed. The thermal efficiency is the fraction 
of the capacitive discharge energy which is channeled into heating of the gas (i.e., increasing the translational-rotational 
temperature). A Sparklet is attached to the end of a free swinging pendulum. The displacement of the pendulum vs 
time due to a single Sparklet discharge is measured using a laser displacement sensor. A theoretical model for the 
Sparklet thermal efficiency is developed assuming a perfect gas. The thermal efficiency of the Sparklet is calculated by 
fitting the predicted displacement vs time with the experiment. The results indicate less than 10% efficiency depending 
on the capacitance size. The thermal efficiency decreases with larger energy deposition magnitudes. 
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